The aim of this study was using Bayesian approach via Gibbs Sampling (GS) for estimating genetic parameters of production, reproduction and health traits in Iranian Holstein cows. Data consisted of 320666 first- lactation records of Holstein cows from 7696 sires and 260302 dams collected by the animal breeding center of Iran from year 1991 to 2010. (Co) variance components were estimated using a multi-trait animal model analyzed via Gibbs Sampling. After convergence, the highest posterior density region of heritability for milk (MY305), fat (FY305), protein (PY305), age at first calving (AFC), calving interval (CI) and somatic cell score (SCS) were 0.255-0.275, 0.195-0.215, 0.195-0.225, 0.260-0.275, 0.065-0.080 and 0.055-0.075, respectively. Genetic correlations ranged from -0.121 (between FY305 and AFC) to 0.914 (between MY305 and PY305) and for phonotypic correlations, it was from -0.083 (between MY305 and SCS) to 0.929 (between MY305 and PY305. The result of this study showed that production traits and AFC have enough genetic variation to develop breeding programs. The estimated genetic correlations suggest that milk production traits and CI would be affected if increasing milk production is the selection goal. The high genetic correlation between CI with SCS suggests that increasing calving interval trait result in an increased SCS.